Ferromagnetism of cold fermions loaded into a decorated square lattice 
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We investigate two-component ultracold fermions loaded into a decorated square lattice, which 
are described by the Hubbard model with repulsive interactions and nearest neighbor hoppings. 
By combining the real-space dynamical mean-field theory with the numerical renormalization group 
method, we discuss how a ferromagnetically ordered ground state in the weak coupling regime, which 
originates from the existence of a dispersionless band, is adiabatically connected to a Heisenberg 
ferrimagnetic state in the strong coupling limit. The effects of level splitting and hopping imbalance 
are also addressed. 

PACS numbers: Valid PACS appear here 



C/2 

cr; 



I 

o 



> 
o 

o 
o 



X 



I. INTRODUCTION 

Ultracold atoms have spurred intensive research ac- 
tivities since the successful realization of Bose-Einstein 
condensation [Il-Q- Furthermore, by loading a periodic 
potential in the system, the optical lattice system 'g'I is 
realized, which opens a new research arena to study quan- 
tum many-body problems in condensed matter physics 
0, Q. This system is formed by the counterpropagat- 
ing lasers, and thereby has an advantage in controlling 
lattice structure, dimensions, interactions, etc. Recently, 
great progress in reaching low temperatures for optical 
lattices hosting fermions has been achieved, which allows 
for the imaging of Fermi surfaces in three dimensions @ , 
the observation of a Mott insulator [HI , the control of 
superexchange interactions [H, [3, etc. These provide 
an important step toward experimental observations of 
a magnetic phase transition in optical lattices. Theo- 
retically, it has already been discussed how magnetically 
ordered states can be realized in some bosonic |14| - [l9l | 
and fermionic |20l- |28| optical lattice systems. (See also, 
continuous systems |29l432| .) 

In this connection, an optical lattice with a decorated 
square lattice structure is particularly interesting (see 
Fig. Ha)), where a wide variety of physics regarding 
magnetism and superfluidity is included. If A (B and 
C) sites in the system are regarded as copper (oxygen) 
ions, it may be reduced to a two-band model related to 
the high cuprates, where the essence of the d-wave 
superconductivity can be explored. Another interesting 
point related to magnetism is the ground-state property 
at half filling, which originates from the characteristic 
band structure (Fig. [Hb)). Due to the existence of a 
dispersionless band (flat band), between two dispersive 
bands, an infinitesimal repulsive interaction drives the 
system to a ferromagnetically ordered state at half fill- 
ing, as guaranteed by the Lieb theorem [33!|. This theo- 
rem ensures the conservation of the total magnetization 
in the ground state of an electron system on a bipartite 
lattice when the interaction is changed from infinitesi- 



mally small interaction to the strong coupling regime. 
This means that there is no level crossing between the 
ground state and low-lying eigenstates under the evolu- 
tion process, which guarantees that there is no quantum 
phase transition at any repulsive interaction. According 
to this theorem, the decorated square lattice considered 
here has a finite magnetization at any repulsive interac- 
tion. Therefore, this specific optical lattice may provide 
a unique avenue to realize the ferromagnetism in cold 
fermions. However, it is not clear how stable the mag- 
netically ordered ground state is against strength of in- 
teractions, level splitting, hopping imbalance, etc, which 
are the accessible control parameters to drive the ferro- 
magnetic state in an optical lattice system. This also 
provides us with an opportunity to study a fundamen- 
tal and important issue in quantum many-body systems: 
what kind of crossover occurs in such a ferromagnetic 
state, when particle correlations are altered systemati- 
cally from a weak to strong coupling regime. 

To make this point clear, in this paper, we study the 
magnetic properties of two-component ultracold fermions 
on the decorated square lattice at half filling. By 
combining the real-space dynamical mean-field theory 
(R-DMFT) with the numerical renormalization group 
(NRG), we show how the magnetically ordered ground 
state is realized in the system. In particular, we eluci- 
date that there is a smooth crossover from a flat-band 
ferromagnetic ground state in the weak coupling regime 
to a ferrimagnetic ground state in the strong coupling 
Heisenberg limit. We note here that flat band ferromag- 
netism is caused by a dispersionless state located just at 
the Fermi level, by lifting its high spin degeneracy by 
an infinitesimal repulsive interaction. In this sense, it is 
governed by the electrons only around the Fermi level 
(weak-coupling). On the other hand, the Heisenberg fer- 
rimagnetism in the strong coupling limit is stabilized by 
the antiferromagnetic exchange interaction induced via 
a virtual hopping process between the nearest neighbor 
sites, where electrons not only around the Fermi level 
but also in the high energy region contribute to the for- 
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mation of magnetic moments. Therefore, when the in- 
teraction increases, the system shows a crossover from 
weak- to strong-coupUng regimes, which is referred to as a 
crossover between flat-band ferromagnetism and Heisen- 
berg ferrimagnetism in this paper. To be more specific, 
we can say that for infinitesimal strength of interaction, 
only the B and C sublattices are polarized ferromagneti- 
cally with a vanishing magnetization of the A sublattice. 
As the interaction increases, a magnetic moment is in- 
duced gradually on the A sublattice in an opposite direc- 
tion to the ferromagnetic order (ferrimagnetism) where 
the total magnetization of the A, B and C sublattices 
is always conserved. In the following, we will explicitly 
show how this type of crossover behavior emerges via the 
calculation of static and dynamical quantities. 

This paper is organized as follows. We introduce the 
model Hamiltonian and summarize numerical techniques 
in Sec. [TTl The stability of the ferromagnetically ordered 
state is discussed in Sec. IIIII We also address the ef- 
fects of level difference and hopping imbalance. A brief 
summary is given in Sec. IIVI 



II. HAMILTONIAN AND METHOD 

We study ultracold two-component fermions on the 
decorated square lattice shown in Fig. [TJ which are de- 
scribed by the following Hubbard Hamiltonian, 



H 



where cj^icia) creates (annihilates) a fermion with spin a 
at site i, and riia = cl^Cia- In our system, spin "up" and 
"down" just are labels to distinguish the two component 
fermions. Here, tij[~ t-^^{t-^'^)] is the nearest- neighbor 
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FIG. 1: (Color online), (a) Decorated square lattice. Its unit 
cell is shown as the dashed line. The square lattice (A) is 
decorated by the B and C sublattices. (b) Density of states 
for the non-interacting case. Note that there is a (5-function 
peak at tj = 0, which reflects a dispersionless band caused by 
a specific geometry of the lattice. 
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hopping in the x (y) direction, t/ is a repulsive interac- 
tion, and Ai is the on-site energy at the ith site, which 
depends on the sublattice indices. 

To discuss the ground state properties of the system, 
we use the R-DMFT, which is regarded as an extended 
version of dynamical mean- field theory (33 - [36j . In this 
method, the lattice problem is mapped onto a set of single 
impurity problems embedded in a self-consistent bath, 
where local particle correlations can be taken into ac- 
count precisely. This method [13] has successfully been 
applied to some correlated systems such as a surface 
an interface of Mott insulator and band insulator 
and cold fermions with a confining potential |2l], l40l 
In general, intersite correlations should be crucial to 
discuss the paramagnetic Mott insulating state in low- 
dimensional systems, which cannot be fully taken into 
account in the R-DMFT method. However, according to 
the Licb theorem, the ground state of the Hamiltonian 
eq. (HI) is always ordered except for limiting cases U — 
or R = 0, where R — t^^ /t^^ . Therefore, it is expected 
that the R-DMFT, which incorporates the sublattice de- 
pendence of particle correlations, can correctly describe 
the ground-state magnetic properties in this system, at 
least, qualitatively. 

In this paper, we consider a decorated square lattice 
system, where the unit cell is composed of three sites, 
as shown in Fig. 1 (a). If we exploit the naive DMFT 
method, in which intersite correlations are completely 
neglected, it is difficult to describe magnetically ordered 
states. However, by using the R-DMFT, which maps 
the original lattice to three kinds of effective impurity 
models in our case, we are still able to properly describe 
the magnetically ordered states. In particular, accord- 
ing to the Lieb theorem mentioned above, our system is 
ensured to have a finite magnetic order at any interac- 
tion strength, so that the R-DMFT provides an efficient 
way to investigate static and dynamical properties of the 
present system. Then the Green function of the effective 
bath Qaa is given by the Dyson equation, 
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-Ai-A-t(k)-E,(w) 



Sqo-(w) 



(2) 



where /i is the chemical potential, Eq is the self-energy 
for the ath sublattice. In the R-DMFT framework, we 
neglect the intersite correlations, which means that S is 
given by a diagonal matrix. The Fourier component of 
hopping t is given by the following 3x3 matrix, 












where k is the wave vector in the reciprocal lattice space 
and a{= 1) is lattice spacing. We note that the self- 
energy does not depend on fc, but on the sublattice and 
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spin indices. This allows us to deal with intersite cor- 
relations to a certain extent through the site-dependent 
self-energies. The effective bath for each sublattice is ex- 
plicitly given as, 



E 



CB.a(^)Cc,a(^ ) 

det{k, w) 



H -1 



+ i;a,<t(w) 



(3) 



E 



det(k, w) 



CA,a(a^)CB,aH-eABW 

det(k, w) 



+ SB,cr(l^) 

(4) 
+ Sc,o-(w) 
(5) 



and 



det(k, uj) 



Ca,£t(w)Cb,ct(w)Cc,<t(w) 

-eic(k)CB,.(c.)-4B(k)Cc,.(^) (6) 



where eAB(k) = -2t^^ cos(fcj^), eAc(k) = -2^^*^ cos(fc:r) 
and Ca,cr(w) = u) + n — Aa — Sqo-(w). By solving the 
three impurity problems, we obtain the corresponding 
self-energies. Green functions for the effective baths are 
then determined self-consistently through the R-DMFT 
equations ©-([SI). 

To discuss the ground state properties of the model 
system, it is necessary to solve the effective quantum im- 
purity models in a proper way. Note here that the density 
of states in our system has a singular property due to the 
flat band at the Fermi level, so that naive perturbation 
theory in U may fail for the half-filled case. Exact di- 
agonalization and quantum Monte Carlo simulations are 
known to provide numerically exact results for a given 
system size, which are expected to reproduce some of 
the present results. However, these methods may not be 
able to resolve extremely small energy scales properly: it 
is not easy to correctly describe the physical properties 
arising from a delta-function like DOS for a flat band. 
In this respect, we think that Wilson's NRG (explaned 
below) is the most reliable method, which can treat low- 
energy properties with extremely high accuracy for any 
interaction strength, allowing us to cope with a delta- 
function like DOS. Since the other methods are not ac- 
curate enough to check our results at present, we confirm 
the validity of our results by checking them in several 
limiting cases. 

We here use the NRG [i^] as an impurity solver. In 
the NRG, one discretizes the effective bath on a loga- 
rithmic mesh by introducing a discretization parameter 
A(> 1). The resulting discrete system can be mapped to 
a semi-infinite chain with exponentially decreasing cou- 
plings which allow us to discuss the properties in- 
volving exponentially small energy scales quantitatively. 
The resulting chain Hamiltonian is then diagonalized it- 
eratively. An important property of this chain Hamil- 



tonian is, that the hopping between successive sites de- 
creases exponentially with increasing chain length. Since 
on the other hand the number of eigenstates increases 
exponentially with each iterative step, it is impossible to 
keep all the eigenstates computed. However, due to the 
decreasing coupling of the chain sites, only a restricted 
number of low-lying states (= L), usually L = 500-1000 
states, need to be retained at each step. The price one 
pays is the loss of information on high energy states at a 
given NRG iteration, which however can be supplied from 
previous iterations. As already demonstrated by a num- 
ber of previous studies, the NRG is especially efficient 
to deal with low-energy properties with high accuracy, 
which is indeed the case for the present model with a fiat- 
band structure. To ensure that the sum rules for dynam- 
ical quantities are fulfilled, we use the complete-basis-set 
algorithm proposed recently [H, In the NRG cal- 
culations, we use a discretization parameter A = 2 and 
keep L — 800 states at each step. We first perform a 
NRG calculation for solving our impurity problems with 
the effective baths expressed in Eq. (3)-(5), from which 
we obtain the self-energy for each sublattice. The re- 
sulting self-energies are then used, in the next step, as 
input functions to calculate effective baths again. In the 
R-DMFT algorithm supplemented by the NRG, these it- 
erative steps are continued until sufficient convergence is 
achieved for physical quantities. 

In the following, we calculate the particle density 
and the magnetization = — riai) to discuss 

how the magnetically ordered ground state is affected by 
repulsive interactions, level splitting, and hopping imbal- 
ance. We also show the local density of states at each sub- 
lattice. In recent experimental developments, the density 
of states for gas systems becomes an observable quantity 
[is! ]. We think that this technique will be extended to 
fermionic optical lattice systems. This fact encourages 
us, therefore, to expect that local quantities shown here 
will be observed in future experiments, too. For simplic- 
ity, we set the magnetization for the A sublattice to be 
positive {ttia > 0). 



III. FERROMAGNETISM AND CROSSOVER 
FROM WEAK TO STRONG COUPLING 

First, we discuss the nature of the ferromagnetically 
ordered ground state for the model with isotropic hop- 
ping tij = t and zero level splitting = 0. The spon- 
taneous magnetization computed for each sublattice is 
shown in Fig. [51 When the interaction U is introduced, 
the magnetization for the A sublattice is induced gradu- 
ally, while those for the B and C sublattices increase dis- 
continuously from zero to m = —0.25 with infinitesimal 
interaction. This singularity is caused by the flat band 
in the LDOS at the Fermi level (Figs. [Sla) and (f)). 
The resulting ferromagnetism is called flat-band ferro- 
magnetism [33|. As the interaction further increases, the 
magnetization for each sublattice increases gradually and 
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FIG. 2: (Color online). Spontaneous magnetization for each 
sublattice: squares (triangles) correspond to the A (B and C) 
sublattice. Crosses represent the total magnetization per unit 
cell, rritot ~ —0.5, which does not depend on the interaction, 
in accordance with the Lieb theorem. 



approaches the saturated values in the strong coupling 
regime, where the fermions are completely localized and 
the system should be described by the antiferromagnetic 
Heisenberg model having the nearest neighbor exchange 
coupling, J ~ P/U. In the decorated square lattice we 
consider here, an additional site is added at the center 
of each bond of the square lattice. Thus, the unit cell 
is composed of three sites (as shown in Fig. 1(a) ), so 
that in the strong coupling limit the system should be in 
a ferrimagnetic state. From a slightly different point of 
view, we can also regard our system as a 1/4 depleted 
square lattice, which naturally explain the emergence of 
a partially depleted antiferromagnetism, i.e. ferrimag- 
netism. Therefore, even though our lattice is bipartite, 
the ferrimagnetic ground state, instead of the antiferro- 
magnetic state, is realized with staggered magnetizations 
(niA ^ 0.5 and mB,c ^ —0.5) in the strong coupling 
Heisenberg limit. In this way, there is a smooth crossover 
from the flat-band ferromagnetism to the Heisenberg- 
type ferrimagnetism when the strength of U is varied. 
We have confirmed that the total magnetization of the 
unit cell is always preserved, rritot = —0.5, during the 
change of the interaction strength, which is consistent 
with the Lieb theorem. 

We can clearly see in the LDOS (Fig. [3]) how such a 
crossover occurs from the weak to strong coupling regime. 
As mentioned above, at U = (Figs. [3] (a) and (f)), 
the LDOS has the flat band in the B, C sublattices at 
the Fermi level, but not in the A sublattice. When an 
infinitesimal interaction is introduced, the splitting of the 
flat band is induced at the Fermi level (Fig. |3](g)), which 
causes the sudden increase in magnetizations for the B 
and C sublattices, in contrast to the smooth increase for 
the A sublattice (Fig. [3] (b)). The weight of the flat 
band is a third of the total weight in the DOS, resulting 
in the total magnetization rritot = —0.5 per unit cell in 
accordance with the computed results in Fig. [5] In this 
way, the ferromagnetism in the weak coupling regime is 
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FIG. 3: (Color online). Local density of states for the A 
sublattice (left panel) and the B, C sublattices (right panel): 
U/t = 0.0, 1.4, 2.8, 5.6, and 11.2 (from top to bottom). 
Solid (dashed) lines represent the LDOS for spin-up (-down) 
fermions. 



dominated by the fermions around the Fermi level in the 
B and C sublattices. As U increases further, the flat-band 
structure is gradually smeared for the B and C sublattices 
(Figs. El (h) and (i)), and the LDOS for A, B and C 
sublattices start to feature simflar profiles (Figs. [3] (d) 
and (i)). Eventually, in the strong coupling regime, where 
the magnetization for each sublattice is almost saturated, 
all the fermions in the whole energy region contribute to 
the formation of the Heisenberg-type ferrimagnetic order. 
Correspondingly, the LDOS shows similar shapes for the 
A, B and C sublattices (Figs. |3](e) and (j)). 

We next consider the effect of the level splitting be- 
tween the A and (B, C) sublattices, A(= Aa — Ab,c), in 
the half-filled system. Such a splitting may be manipu- 
lated experimentally by introducing a site-dependent po- 
tential. By performing the R-DMFT calculations with 
the total number of particles fixed, we obtain the results 
for a specific choice of the interaction, U/t — 5.6, which 
are shown in Fig. S) Note that the obtained magnetiza- 



5 



m 




FIG. 4: (Color online), (a) Spontaneous magnetization 
and (b) particle density of each sublattice when U/t = 5.6. 
Squares (triangles) represent the magnetization for the A (B 
or C) sublattice. Crosses represent the total magnetization 
per unit cell. 



tions are symmetric about A = 0. When A = 0, the sys- 
tem is in the region of ferrimagnetic ground state, where 
the sublattice magnetization is large enough 0.3 
and \mB,c ~ 0.4|). The introduction of the level split- 
ting A decreases (increases) the particle density for the 
sublattice A (B and C), as shown in Fig. [5l At the same 
time, the magnetization for each sublattice is monotoni- 
cally decreased, and the effect of onsite correlations be- 
comes less important. When |A|/t = 11.2, the sublattice 
magnetization for the A sublattice almost vanishes, but 
is finite for the B and C sublattices. This implies that 
the crossover occurs from the ferrimagnetic state to the 
fiat-band ferromagnetic state when |A| is increased. The 
crossover is indeed confirmed by observing LDOS cal- 
culated for the same parameters shown in Fig. [5] It 
is seen that a J-function like peak characteristic of the 
flat-band ferromagnetism appears in the vicinity of the 
Fermi level only in the LDOS for the B and C sublattices 
at |A|/t ^ 11.2. Namely, the magnetization is governed 
by fermion states near the Fermi level for the B and C 
sublattices. 

We would like to comment here on a technical prob- 
lem in our numerical calculations, which arises when we 
numerically fix the total number of particles in the pres- 
ence of large level difference A. When the level difference 
is introduced, the DOS does not preserve particle-hole 
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FIG. 5: (Color online). Local density of states for the A 
sublattice (left panel) and the B, C sublattices (right panel) 
when U/t = 5.6: A/t =-11.2, -5.6, 0.0, 5.6, and 11.2 (from 
top to bottom). Solid (dashed) lines represent the LDOS for 
up (down) spin fermions. At A/t = 11.2 (-11.2), a 5-function 
like peak appears for up-spin (down-spin) fermions in the B,C 
sublattices. 



symmetry, which was used efficiently to fix the particle 
number at half-filling for A = 0. Therefore, the introduc- 
tion of finite A requires more elaborate numerical calcu- 
lations to keep the filling factor unchanged with sufiicient 
accuracy. Thus DMFT calculations typically take much 
longer time to adjust the chemical potential properly in 
each iteration step. For larger level splitting, the ground 
state is very sensitive to the change in the chemical po- 
tential, so that this sort of problem becomes more serious. 
Nevertheless, the half-filled system we address here satis- 
fies the condition of the Lieb theorem for any value of A, 
which ensures that the magnetic state should be stable 
even for very large |A|, although the ground state may 
become very sensitive to a small amount of particle or 
hole doping. The results deduced numerically for large 
|A| are indeed consistent with the arguments based on 
the Lieb theorem. 

To discuss the crossover between two extreme limits 
of the magnetically ordered state, we systematically per- 
form similar calculations for various choices of the level 
splitting and the interaction. The obtained results for 
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FIG. 6: (Color online). Spontaneous magnetization for (a) 
A sublattice (b) B and C sublattices as functions of U/t and 
A/t. Bright (dark) region has a large (small) magnetization. 
It is seen that a crossover appears around the lines of [/ ~ | A| 
and U ~ W. 



the magnetization are shown (shading plots) in Fig. [51 
When U/t is large and |A|/t is small, the onsite Coulomb 
interaction plays an important role in stabilizing the fer- 
rimagnetically ordered ground state with \ma\ ~ 0.5, as 
mentioned above. On the other hand, when the level 
splitting is large, the magnetization for the A sublattice 
vanishes and those for the B and C sublattices are close 
to half of full polarization, mB,c ^ —0.25, characteristic 
of flat-band ferromagnetism. It is seen in Fig. [5] that 
the crossover between two types of magnetically ordered 
states appears around |A//7| ~ 1 and U/W ~ 1, where 
W{= AV2t) means the bandwidth at A = 0. Let us com- 
ment on the crossover as a function of the level difference 
A found here. This crossover is related to the original one 
as a function of U discussed above. In fact, the physics 
emerging in both cases are essentially the same. The cru- 
cial point is that a flat band always exists for any value 
of the level difference A at [/ = 0, and it consists of a 
specific combination of the states belonging to the B and 
C sublattices. For finite A, the strength of the effective 
interaction should be compared with the value of A: for 
a given value of U, larger (smaller) A effectively drives 
the system to a weak (strong) coupling regime. There- 
fore the change in the level difference causes the change 
in the effective interaction, which in turn gives rise to 
a crossover from weak- to strong-coupling regimes. One 
can indeed see the formation of a delta-function like peak 
in Fig. 5(i) and (j) for large A, which is characteristic 




FIG. 7: (Color online). Spontaneous magnetizations as func- 
tions of Ri= t^'^/t^^) when U/t^^ = 2.8. 



of the weak-coupling regime (fiat-band ferromagnetism). 
We note that the Lieb theorem even in this case guaran- 
tees no transition under the change of A 

We finally discuss the effects of hopping imbalance 
R{— t^'^/t^^), which may be tunable by laser tech- 
niques. When R is decreased from unity, hopping of 
fermions becomes spatially anisotropic. For example, in 
the small R region, the system can be regarded as the 
one-dimensional Hubbard chains in y-direction that are 
weakly coupled to isolated sites (C sublattice). In the 
case, the absolute values of the magnetizations for the 
A and B sublattices seem to approach the same finite 
value (Im^^sl ^ 0.15) as R decreases, as shown in Fig. 
[T] On the other hand, the spins for the C sublattice 
become fully polarized, since those spins become almost 
free. Note that as far as i? 7^ 0, the system satisfies 
the condition of the Licb theorem, and thereby the total 
magnetization, mtot — —0.5, should not be altered during 
the change of i?, as is indeed the case in Fig. [T] In the 
special case of i? = 0, the system is completely decou- 
pled into two subsystems consisting of Hubbard chains 
and isolated sites. The Hubbard chain realized at i? = 0, 
which consists of A and B sublattices, corresponds to a 
ID system with periodic boundary conditions. However, 
in the DMFT the information of the system is incorpo- 
rated through the shape of the density of states. So, it is 
rather difficult to deduce physical properties correctly in 
this limiting case. This point thus remains an open issue 
for future studies. In the Hubbard chain, a magnetic long 
range order does not appear even at absolute zero, and 
the associated elementary spin excitations are gapless. 
Therefore, it might be expected that on the introduction 
of infinitesimal R, isolated free spins (C sublattice) affect 
the Hubbard chains to induce finite sublattice magneti- 
zations in the A and B sublattices. Since our mean-field 
calculations may not be appropriate to discuss the weakly 
coupled chains, magnetic properties around ~ should 
be examined by incorporating intersite correlation more 
precisely, which will be addressed elsewhere. 

Before closing this section, we would like to comment 
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on the validity of our study. In the framework of the 
R-DMFT, we have fuhy incorporated onsite correlations, 
but have not dealt with intersite correlations sufficiently. 
This indeed leads to the sublattice magnetizations that 
are somewhat overestimated, particularly in the strong 
coupling regime: the; magnetization for each sublattice 
takes almost saturated value, \mA,B.c\ ^ 0.5, which 
should be suppressed if we introduce quantum fluctu- 
ations originating from intersite correlations more pre- 
cisely. However, we have confirmed that our calculations 
always satisfy the condition of the Lieb theorem that 
the total magnetization is constant, mtot = —0.5, in the 
whole parameter space. Note that the R-DMFT calcula- 
tion of the ground state is performed independently for 
various choices of interactions, from which we figure out 
how large the ground-state magnetization is. The above 
confirmation is very important because the R-DMFT is a 
mean field theory and therefore it is not a priori obvious 
whether this method is applicable to the fiat-band system 
considered here. We have also revealed that the crossover 
in the density of state interpolates correct behaviors both 
at weak and strong coupling limits. This result also en- 
courages us to believe the validity of the method. There- 
fore, we believe that our R-DMFT study can describe the 
essential properties of the system at least qualitatively 
in the whole parameter region, although intersite mag- 
netic correlations should be taken into account in order 
to discuss the results more quantitatively in the strong 
coupling regime and for i? — ^ 0. 

IV. SUMMARY AND DISCUSSIONS 

We have investigated the nature of the ferromagnet- 
ically ordered ground state in the half-filled Hubbard 
model on a decorated square lattice by means of R- 
DMFT and NRG. By varying the interaction strength, 
the level splitting, and the hopping imbalance, we have 
found that the ferromagnetic state, which originates from 
a fiat-band structure, is adiabatically connected to the 
Heisenberg ferrimagnetic state in the strong coupling 
regime. The magnetic ground state becomes more sta- 
ble as U increases from zero, but beyond the crossover 
regime it is governed by the Heisenberg exchange inter- 
action, ^ /U , whose magnitude decreases as U further 
increases. Therefore, the energy gain forming the mag- 
netically ordered state is largest in the crossover regime, 
which implies that the magnetic ground state should be 



most stable in the crossover regime. To realize the above 
magnetic state experimentally in optical lattice systems, 
it is important to examine how stable it is at finite tem- 
peratures. In the pure 2D case discussed here, a finite- 
temperature phase transition does not occur according to 
the Mermin- Wagner theorem. If we consider a quasi-2D 
system (e.g. a collection of weakly coupled layers, etc), 
however, a phase transition can occur at a finite tem- 
perature. The corresponding transition temperature is 
expected to be highest in the crossover regime, according 
to the above discussions on the stability of the ground 
state. This fact encourages experimental investigations 
on a quasi-2D optical lattice system, for which the system 
parameters are tuned in the crossover regime. 

Another point to be noticed in optical lattice systems 
is that there is a confining trap for ultracold atoms, which 
superposes a smoothly varying harmonic potential to the 
periodic potential for the decorated square lattice. In 
a local density approximation, it can be regarded as a 
spatially modulated chemical potential. Therefore, for 
small U, the area which satisfies the half-filling condition 
relevant to forming a local Mott state, if it exists, may 
be very small around the center of the system. As U 
increases, the area of local Mott state is extended, and 
at the same time, the ferromagnetic state becomes more 
stable, changing its character from a fiat-band type to a 
Heisenberg type gradually. Therefore, it is expected that 
for reasonably large interactions around the crossover 
regime, the ferromagnetically ordered state may be ob- 
served in the optical lattice. Fortunately, this conclu- 
sion is consistent with the temperature effects mentioned 
above. Detailed discussions on a trap potential will be 
addressed elsewhere. 
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